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INTRODUCTION 

Humidity variability at the top of the marine boundary layer (MBL) and 
in the free troposphere has been examined using a variety of measurements 
taken on and around San Nicolas Island (SNI) during the FIRE IFO in July, 
1987. Doppler wind profiler reflectivity recorded at two minute time 
resolution has provided the most continuous record and detail of small scale 
humidity fluctuations. Rawinsonde data was available from both an island 
site and the research vessel Point Sur . The information extractable from 
these sources is somewhat limited due to the frequency of launches (3-4/day 
at SNI and 6/day on the Point Sur). Some additional data was available from 
instrumented aircraft although scheduling flights in the neighborhood of the 
island was difficult due to restrictions on the air space. Other relevant 
data was collected at SNI near the radar and rawinsonde launch sites. A 
continuous record of cloud base altitude was logged by a ceilometer. Doppler 
acoustic sounder (sodar) reflectivity data provided a good record of 
inversion height. The sodar also monitored turbulent temperature 
fluctuations in the MBL. A small ground station recorded hourly averages of 
solar irradiance and downward longwave irradiance. 

This paper describes analysis in progress of the various data sets 
described above for two adjacent two day periods from 11 July to 14 July. 

The earlier period was chosen because the marine inversion was unusually high 
and there was increased frequency of rawinsonde launches at SNI . The later 
period was chosen because of the significant descent with time of an elevated 
inversion indicated by the radar data. Throughout the four day period, but 
especially in the first half, the turbulent humidity structure calculated 
from Doppler radar reflectivity shows excellent agreement with humidity 
profiles evaluated from rawinsonde data. 

INTERPRETATION OF DOPPLER RADAR REFLECTIVITY 

Penn State's Doppler wind profiler used at SNI operates at a frequency 
of 404.37 MHz. For FIRE, the pulse width was set at 1 ^s to provide the best 
possible range resolution while satisfying the receiver bandwidth constraint. 
The pulse repetition frequency was 10,000 Hz. Processing of radar 
reflectivity is accomplished in the following manner. The complex video of 
288 consecutive pulses are coherently averaged (time domain integration). A 
set of 64 of these integrations is processed by a fast fourier transform 
(FFT) to produce one spectrum. Incoherent averaging is then performed on 8 
of these variance spectra to produce the two minute averaged spectrum from 
which the signal/noise ratio (SNR) is extracted. 

Since radar backscatter occurs in a volume (V ) defined by the pulse and 
beam widths, volume reflectivity (r?) replaces the more familiar backscatter 
cross -section (a), where r}=o /V . The volume reflectivity is related to the 
SNR by the spectral radar equation (VanZandt et al . , 1978) 
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where c is the speed of light, k is Boltzmann's constant, B is the bandwidth 
of the integrating filter, a is the combined antenna/line efficiency, T 
a nd T are noise temperatures for cosmic noise and the receiver, P is°the 
transmitted power, is the pulse repetition frequency, A is the antenna 
area, x is the off vertical beam axis angle, R is the rangE, and AR is the 
range resolution. Based on comparison with a previously calibrated 50 MHz 
Doppler wind profiler located in central Pennsylvania, the overall antenna 
efficiency was chosen to be 0.18. Additional calibrations of the antenna 
using aircraft data are currently underway. At A-0.742 m, qT «T so to a 
first approximation cosmic noise interference may be neglecteS. ftith the 
appropriate values inserted, (1) becomes 


r? - 7.67x10' 28 R 2 SNR. 


( 2 ) 


The backscatter intensity of the radar's signal depends on the mean 
refractive index gradient existing in the scattering volume. Therefore, the 
refractive index structure function parameter (C ) can be written in terms 
of the radar's volume reflectivity (Ottersten, 1?69) 

C n 2 - (r,/0.38)A (1/3) . (3) 

Insertion of A-0.742 m and (2) into (3) yields 

C n 2 - 1 . 8 3x10 2 7 R 2 SNR. ( 4 ) 


Both temperature and humidity fluctuations may contribute to variations 
i^the refractive index within the scattering volume. In radar studies of 
at l° w altitudes, the contribution due to temperature fluctuations is 
usually ignored (Wesley, 1976). Conversion from C to C 2 is then 
(following Burk, 1980) n 9 

Cq 2 - C n 2 (1667T 2 /P) 2 , (5 ) 

where T is the temperature and P is the pressure (in mb) within the 
scattering volume. Combinbin| (5) and (4) yields the final form of the 
equation used to determine from Doppler radar reflectivity; 

C^ 2 = 5.08x10 " 21 R 2 SNR(T 2 /P) 2 . (6) 

RAWINSONDE DATA 


During the FIRE IFO, Colorado State University operated a cross-chain 
Loran atmospheric sounding system (CLASS) at SNI . Scientists from the Naval 
Postgraduate School launched VIZ Loran- type sondes from the Point Sur . For 
additional information on the CLASS system and details on processing of the 
SNI rawinsonde data, see Schubert, et al . (1987a). Water vapor mixing ratio 
(q) was computed in the standard way. 
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RESULTS 


Fig. 1 shows a time-height cross-section of log(C ) contours (g/kg) 
from (6) for July 11-12. In this and the following figures, time in ^ours is 
measured continuously starting at Oz on 11 July. Higher values of C 
correspond to regions of increased scattering from turbulent humidit^ 
fluctuations . ^ In regions of active turbulence, the structure function 
parameter (C ) for a particular variable (x) is proportional to the square 
of the vertical gradient of x (Fairall et al , 1988) 

c x 2 - 1.6 7x (C u 2 /N 2 )(ax/az) 2 , (7) 

2 

where x could be T, q, or n; C , a measure of local turbulence, is the 

velocity structure function parameter, N is the Brunt-Vaisala frequency; and 

7 DO . 3 is a constant. Thus, C should correlate well with vertical humidity 
x q J 

gradients . 

For comparison, Fig. 2 shows q contours (g/kg) derived from rawinsonde 
data for the same time period shown in Fig. 1. The swath of high Jq/Jz 
across the middle of Fig. 2 is evidence of the MBL inversion height 
decreasing from about 1000 m at the beginning of 11 July to about 730 m at 
the end^of 12 July. Note that the humidity gradients argj well correlated 
with C except for a period from 23z to 29z in which C reaches a local 
minimuS while dq/dz approaches its maximum value. ^ 

This apparent disagreement with the relationship in (7) can be analyzed 
using additional data obtained from the remote sensors located at SNI . A 
record of cloud base measured by the ceilometer (Schubert et al . , 1987b), is 
shown in Fig. 3. Clear sky is depicted by a point at 960 m. The period 
labeled 23z to 27z shows a sharp decrease in cloud base, suggesting a period 
of strong subsidence which eventually leads to clearing just after 27z. Wind 
speed contours from merged sodar/wind profiler data sets (Syrett, 1988) are 
shown in Fig. 4. Here, the period from 23z to 29z is characterized by light 
winds and negligible vertical wind shear. Thus for the period in question, 
one may conclude that although a strong humidity gradient exists, the 
mechanisms that would normally generate turbulence in the region of the 
gradient such as cloud top entrainment a^d vertical wind shear, are not 
present, resulting in lower values of . 
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Fig. 1 Time-height cross-section of log(C ) contours (g/kg) calculated from 
Doppler radar reflectivity. ^ 

Fig. 2 Time-height cross-section of specific humidity contours (g/kg) from 
rawinsondes launched at SNI. Data was compiled by Schubert et al . , (1987a). 

Fig. 3 Filtered cloud base versus time from ceilometer data collected at 
SNI during FIRE (after Schubert et al . , 1987b). Clear sky is depicted by a 
point at 960 m. 

Fig. 4 Time-height cross-section of wind speed contours (m/s) derived from 
merged sodar/wind profiler data sets (after Syrett, 1988). 
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Fis- 3 Filtered cloud base versus time from ceilometer data collectad at 
SHI during FIRE (aftar Schubert at al., 1987b). Claar sky is depicted by a 
point a 960 m. 
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Fig. 4 Time-height cross-section of wind spaed contours (m/a) derived from 
merged sodar/wind profiler data sats (after Syrett, 1988). 
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Fig. 1 Tln»-h»ight cro»*-**ctlon of log(C ) contour* (g/kg) calculated froa 
Doppler radar reflectivity. 



TIME (G.M.T.) 


Fig. 2 Time-height cross-section of specific humidity contours (g/kg) from 
rawinsondes launched at SNI . Data was compiled by Schubert at al., (1987a). 
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